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We update the non-perturbative function of tlie Collins-Soper-Sterman (CSS) 
resummation formalism which resums the large logarithmic terms originating 
from multiple soft gluon emission in hadron collisions. Two functional forms 
in impact parameter (6) space are considered, one with a pure Gaussian form 
with two parameters and another with an additional linear term. The re- 
sults for the two parameter fit are found to be gi = 0.24lQQy GeV^, §2 = 
OMto^ol GeV^. The results for the three parameter fit are gi = O.l^tom GeV^, 
c,2 = 0^48t[J;[]^ GeV^, and 53 = -0.58toit GeY'^. We also discuss the potential 
of the full Tevatron Run 1 Z boson data for further testing of the universality 
of the non-perturbative function. 



Abstract 



PACS numbers: 12.15.Ji,12.60.-i,12.60.Cn,13.20.-v,13.35.-r 



1 Introduction 



It is a prediction of the theory of Quantum Chromodynamics (QCD) that at hadron 
colhders the production of Drell-Yan pairs or weak gauge bosons {W^ and Z) will 
generally be accompanied by gluon radiation. Therefore, in order to test QCD theory 
or the electroweak properties of vector bosons, it is necessary to include the effects 
of multiple gluon emission. At the Fermilab Tevatron (a pp collider), we expect 
about 2 X 10^ W"^ and 6 x 10^ Z bosons produced at = 1.8 TeV, per 100 pb"^ 
of luminosity. This large sample of data is useful (i) for QCD studies (with either 
single or multiple scales), (ii) as a tool for precision measurements of the W boson 
mass and width, and (iii) as a probe for new physics (e.g., Z'). Achievement of these 
physics goals requires accurate predictions for the distributions of the rapidity and 
the transverse momentum of W^/Z bosons and of their decay products. 

Consider the production process /ii/i2 — ^ VX. Denote Qt and Q to be the trans- 
verse momentum and the invariant mass of the vector boson V, respectively. When 
Qt ~ Q, there is only one hard scale and a fixed-order perturbation calculation is 
reliable. When Qj- <^ Q, there are two hard scales and the convergence of the con- 
ventional perturbative expansion is impaired. Hence, it is necessary to apply the 
technique of QCD resummation to combine the singular terms in each order of per- 
turbative calculation, which yields: 

+a|(L5 + L^) + 4(L^ + L6) + ... (1) 
+«|(L3 + l2) + 4(L5 + L^) + ... 
+a|(L+ l) + a|(L3 + L2) + ... }, 

where is the strong coupling constant, L denotes \n{Q'^/Q'^) and the explicit coef- 
ficients multiplying the logs are suppressed. 

Resummation of large logarithms yields a Sudakov form factor [|l], Q and cures 
divergences as Qt —>■ 0. This resummation was pioneered by Dokshitzer, D'yakonov 
and Troyan (DDT) who performed an analysis in Q^-space which led to a leading- 
log resummation formalism [Q]. Later, Parisi-Petronzio showed [0] that for large Q 
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+a|(L+ 1) 



the Qt region can be calculated perturbatively by imposing the condition of 
transverse momentum conservation, 

in 6-space {b is the impact parameter, which is the Fourier conjugate of Qt)- Their 
improved formalism also sums some subleading-logs. They showed that as Q ^ oo, 
events at Qt ~ may be obtained asymptotically by the emission of at least two 
gluons whose transverse momenta are not small and add to zero. The intercept at 
Qt = is predicted to be 

da 

dQT Q^^o 

where rjo = Alia 1 + ^ with A = 12Cf/{S3 — 2nf), and r/o — 0.6 for Uf = 4 and 
Cp = 4/3. Collins and Soper extended this work in 6-space and applied the 
properties of the renormalization group invariance to create a formalism that resums 
all the large log terms to all orders in a^. 

Although various formalisms for resumming large In^Q^/Q"^) terms have been 
proposed in the literature [|, ^, we will concentrate in this paper on the formalism 
given by Collins, Soper and Sterman (CSS) 0, which has been applied to studies of 



-0 7^ , (3) 



the production of single [|^ |], |, [T^ and double |[TT| weak gauge bosons as well as 
Higgs bosons |12| at hadron colliders. 



2 Collins-Soper-Sterman Resummation Formalism 

In the CSS resummation formalism, the cross section is written in the form 

where y is the rapidity of the vector boson V, and the parton momentum fractions 
are defined as Xi = e^Q/y^ and X2 = e'^Q/^/S with as the center-of-mass 
(CM) energy of the hadrons hi and /i2- In Eq. (p, y is the regular piece which can 
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be obtained by subtracting the singular terms from the exact fixed-order result. The 
quantity W satisfies a renormalization group equation with the solution of the form 



Here the Sudakov exponent is defined as 



S{Q,b,C,,C2)= I 
Jc 



A («,(/!), Ci) In (^^] +B (as(7l), Ci, C2) 



(6) 



and the xi and X2 dependence of W factorizes as 

Here, Cjh is a convolution of the parton distribution function (fa/h) with calculable 
Wilson coefficient functions {Cja), which are defined through 

CjhiQ, b,x) =J2j^ y Cja 6, ^ = y , fa/h (^^, /i = y ) • (8) 

The sum on the index a is over incoming partons, j denotes the quark fiavors with 
(electroweak) charge ej, and the factorization scale fi is fixed to be C3/6. A few 
comments about this formalism: 

• The A, B and C functions can be calculated order- by-order in a^. 

• A special choice can be made for the renormalization constants Ci so that the 
contributions obtained from the expansion in of the CSS resummed calcu- 
lation agree with those from the fixed-order calculation. This is the canonical 
choice. It has Ci = C3 = 2e~'''^ = 60 and C2 = Ci/bo = 1, where ■je is Euler's 
constant. 

• 6 is integrated from to 00. For b ^ 1/Aqcd, the perturbative calculation is no 
longer reliable. In order to account for non-perturbative contributions from the 
large b region this formalism includes an additional multiplicative factor which 
contains measurable parameters. 

We refer the readers to Ref. [^] for a more detailed discussion on how to apply the 
CSS resummation formalism to the phenomenology of hadron collider physics. 



2.1 The Non-perturbative Function 

As noted in the previous section, it is necessary to include an additional factor, usually 
referred to as the "non-perturbative function", in the CSS resummation formalism in 
order to incorporate some long distance physics not accounted for by the perturbative 
derivation. Collins and Soper postulated 

Wfkib) = W,MW;i^{b) , (9) 

with 

K = I ^ , (10) 

SO that h never exceeds hmax and Wjiip^) can be reliably calculated in perturbation 
theory. (In numerical calculations, hmax is typically set to be of the order of 1 GeV^^.) 
Based upon a renormalization group analysis, they found that the non-perturbative 
function can be generally written as 



W-i {b,Q,Qo,Xi,X2) = exp 



-Fi(6)ln ~ Fj/h,{xi,b) - Fj,/,^^{x2,b) 



, (11) 



where Fi, Fj/h^ and Fj^//^^ must be extracted from data with the constraint that 

Wf/ib = 0) = 1. 

Furthermore, Fi only depends on Q, while Fj/^^ and Fj^//^^ in general depend on xi 
or X2, and their values can depend on the flavor of the initial state partons (j and k 
in this case). Later, in Ref. |T3[, it was shown that the Fi{b) In (^q^^ dependence is 
also suggested by infrared renormalon contributions to the Qt distribution. 

2.2 Testing the Universality of 

The CSS resummation formalism suggests that the non-perturbative function should 
be universal. Its role is analogous to that of the parton distribution function (PDF) 
in any fixed order perturbative calculation, as its value must be determined from 
data. The first attempt to determine such a universal non-perturbative function was 
made by Davies, Webber and Stirling (DWS) [jl4[ in 1985 who fit data available at 



that time to the resummed piece (the Wj]^ term) Duke and Owens parton distribution 
functions . Subsequently, the DWS results were combined with a NLO calculation 
Ill6| by Arnold and Kauffman in 1991 to provide the first complete CSS prediction 
relevant to hadron collider Drell-Yan data. In 1994, Ladinsky and Yuan (LY) [|1^ 
observed that the prediction of the DWS set of deviates from R209 data {p+p — > 
12+ 12- + X at \^ = 62GeV) using the CTEQ2M PDF p|. In order to incorporate 
possible ln(r) dependence, LY postulated a model for the non-perturbative term, 
which was different from that of DWS, as 

Q 



Wj-^ (6, Q, Qo, xi, xs) = exp 



-gib - g2b In 



2Qo 



- gigsbln {100x1X2) 



(12) 



where X1X2 = t. A "two-stage fit" of the R209, CBF-Z (4pb~^ data) and E288 



{p + Cu) data gave ||T7| 



gi = 0.1lt'ofs GeV^ = 0.58l°i GeV^ 



, g, = -1.5i°:i GeV-\ 



for Qo = 1.6 GeV and bmax = 0.5 GeV~^. [] The purpose of the project described here 
is to update these non-perturbative parameters using modern, high-statistics samples 
of Drell-Yan data and to incorporate a fitting technique which will track the full error 
matrix for all fitted parameters. Our results are given in the following sections. 



3 Fitting Procedure 

3.1 Choice of the Parametrization Form 

At the present time, the non-perturbative functions in the CSS resummation formal- 
ism cannot be derived from QCD theory, so a variety of functional forms should be 
studied. The only necessary condition is that W^^{b = 0) = 1. For simplicity, we 
consider only two typical functional forms for Wj^^{b,Q,QQ,Xi,X2) in b space: (i) 
2-parameter pure Gaussian form [DWS form]: 

Q ^ 



exp 



-gi - 92 In 



2go. 



(13) 



FORTRAN coding error in calculating the parton densities of the neutron led to an incorrect 
value for 33. 
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with Qo = 1.6 GeV; and (ii) a 3-parameter form [LY form]: 

expl -5(1 - 5(2 In (-^) - [gigs\n{100xiX2)]b 




(14) 



with a logarithamic x-dependent third term which is hnear in b. This is equivalent 
to ln(r/ro) for Tq = 0.1. 

Both forms assume no flavor dependence for simplicity. In addition to fitting 
for the non-perturbative parameters, gi, g2, andg^, the overall normalizations were 
allowed to float for some fits. One can also study another pure Gaussian form with 
similar x dependence such as 



However, we find that current data are not yet precise enough to clearly separate 
the g2 and g^ parameters within this functional form and so it is not considered 
here. We also tested a few additional functions which did not incorporate additional 
parameters, but did not find any clear advantage to them when fitting the current 
Drell-Yan data. However, as to be shown later, the Run IW/Z data at the Tevatron 
are expected to determine the g2 coefficient with good accuracy, and these data can 
be combined with the low energy Drell-Yan data to further test various scenarios for 
X dependence and ultimately, universality. 



In order to determine the non-perturbative functions discussed above, we need to 
choose experimental data sets for which the contribution to the non-perturbative piece 
dominates the transverse momentum distributions. This suggests using low energy 
fixed target or collider data in which the transverse momentum {Qt) of the Drell- 
Yan pair is much less than its invariant mass (Q). Because the CSS resummation 
formalism better describes data in which the Drell-Yan pairs are produced in the 
central rapidity region (as defined in the center-of-mass frame of the initial state 
hadrons) we shall concentrate on data with those properties. Based upon the above 
criteria we chose to consider data shown in Table |^. We have also examined the E772 




3.2 Choice of the Data Sets 
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Experiment 


reference 


Reaction 


y/S GeV 


(r) 


6N 


R209 


m 


p + p ^ + X 


62 


~ 0.1 


10% 


E605 


m 


p + Cu-^ + X 


38.8 


~ 0.2 


15% 


CDF-Z 


m 


p + p^ z + x 


1800 


~ 0.05 




E288 


m 


P+CU-* /i+yU" + X 


27.4 


~ 0.2 


25% 



Table 1: Drell-Yan data used in this analysis. Here, 6N is the published normalization 
uncertainty for each experiment. The CDF data were from Tevatron collider Run 
of 4 pb-^ 

data [|l^, from the process p + — > /i^/i^ + X at = 56.6 GeV, and found that 
it was not compatible in our fits with the above data, and it is not included in this 
study.0 Except where noted, all of the fits to (71,2,3 were done using the CTEQ3M 
PDF |21[ fits. 

3.3 Primary Fits 

As to be shown later, the E288 data have the smallest errors, and would be expected 
to dominate the result of a global fit. That is indeed the case. When including the 
E288 data in a global fit, we found that the resulting fit required the NORMQ to be 
too large (as compared to the experimental systematic error) for either the E288 or 
the E605 data. Furthermore, the shape of the R209 data cannot be well described by 
the theory prediction based on such a fit. 

3.3.1 Fits ^2,3 

We therefore employed a different strategy for the global fit based on the statisti- 
cal quality of the data. We included the first two mass bins (7 < Q < 8 GeV and 
8 < Q < 9 GeV) of the E605 data and all of the R209 and the CDF-Z boson data, in an 
initial global fit, referred to here as Fit ^2,3. In total, 31 data points were considered. 

We allowed the normalization of the R209 and E605 data to float within their overall 

^ Using the fitted g values to be given below, the theory prediction for the E772 experiment is 
typically a factor of 2 smaller than the data. Similarly, CTEQ fitting of PDF parameters are not 
well fit with these data ||2^ . 

•^Here the quantity NORM is the fitted normalization factor which is applied to the prediction 
curves in all that follow: the data are uncorrected. 
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R209 data (PRL 47 (1981) 12) CTEQ3 
gl = .24, g2= .34. g3= 0., NORM= 1.04 
g1 = .15, g2= .48, g3= -.58, NORM= 1.04 



5<Q<8 GeV 




8<Q<1 1 GeV 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Q, (GeV) 



Figure 1: R209 data, from p + p ^ + X VS = 62 GeV, with an overall 

systematic normalization error of 10%. The curves are the results of Fit A2 and A3 
and are multiplied by the value of NORM, as shown in the figure and described in 
the text. 



CDF Z data { PRL 67 (1991) 2937 ) 

CTEQ3, NORM= 1. 



O g1= .24, g2= ,34, g3= 0, 
• g1 = .1 5, g2= .48, g3= -.58 



10 12. b 15 17.5 20 



Q, (GeV) 



Figure 2: Comparison of 4 pb ^ CDF-Z data at the Tevatron with two different theory 
model predictions. The dots correspond to the results of Fit A2 and A3. 



8 



K605 data (PHD 43 (1991) 2815) 

CTEQ3, s = (38.8 GeV)' 



gl= ,24, g2= .34, g3= 0,, NORM= 0.94 
g1 = ,1 5, g?= .48, g3= -.58, NORM= 0,87 



7<Q<8 GeV 



10.5<Q<1 1.5 GeV 



13,5<Q<18 GeV 



Figure 3: E605 data, from p + Cu — >■ + X at ^/S = 38.8 GeV, with an overall 

systematic normalization error of 15%. The curves are the results of Fit A2 and 
and are multiplied by the value of NORM, as shown in the figure and described in 
the text. 

systematic normalization errors, while fixing the normalization of the CDF-Z data to 
unity. (The point-to-point systematic error of 10% for the E605 data has also been in- 
cluded in the error bars of the data points shown in Fig. ^) In addition to the normal- 
ization factor for each experiment, the fitted parameters of our global fit include the 
coefficients gi^2 and ^'1,2,3 for the 2-parameter and 3-parameter fits, respectively. We 
found that, for Qq = 1.6 GeV and bmax = 1/(2 GeV), both the 2-parameter and the 3- 
parameter forms give good fits, with per degree of freedom about equal to 1.4.Q The 
best fitted central values for Fit A2 are: gi = 0.24 GeV^ , g2 = 034 GeV^ While the 
central values for Fit A3 are gi = 0.15 GeV^ , ^2 = 0.48 GeV^ , gs = -0.58 GeV"\ 
The fitted values for NORM are 1.04 for both Fits A2 and A-,. 



3.3.2 Uncertainties in the Fits 

We have also studied the uncertainties of the fitted g parameters. For the 2-parameter 

fit, the 1(7 error in the plot (with an approximately elliptical contour) gives —0.07 < 

6gi < 0.08 GeV^ and —0.08 < 6g2 < 0.07 GeV^. For the 3-parameter fit, the situation 
^ We scan the values of gi and 52 between and 1, and between —2 and 3. 
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2 PARAMETER FIT 

0.425 
0.4 
0.375 

a. 35 

0.325 
0.3 
0.275 
0.25 

0.15 0.18 0.2 0.22 0.24 0.2S 0.28 0.3 0.32 0.34 

gi 

Figure 4: The error ellipse on the gi and g2 plane from which the errors of the 
2-parameter fit A2 were interpreted. 

is more complicated, as the fitted values of g^s are highly correlated. In order to 
estimate the uncertainties of the fitted g values, we fix the value (at its best fit value) 
of g^s, one at a time, and examine the uncertainties of the other two, in a way similar 
to studying the 2-parameter fit result. We found that 



9i 


fixed: 


-0.05 


< 


Sg2 < 0.04 GeV^ , 






-0.02 


< 


6g-i < 0.26 GeV-^ 


92 


fixed: 


-0.02 


< 


dgi < 0.03 GeV^ , 






-0.08 


< 


5g3 < 0.07 GeV-^ 


93 


fixed: 


-0.02 


< 


6gi < 0.04 GeV^ , 






-0.03 


< 


Sg2 < 0.02 GeV^ ; 



constitute a conservative set of uncertainty ranges. With the understanding of the 
complexity discussed above, we characterize the uncertainties of the fitted g^s in the 
3-parameter form conservatively by their maximal deviations, so that the best fitted 
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3 PARAMETER FIT 



0.51 — 



0.4-9 — 



0.48 — 



O.A-7 — 



O.AS — 



0.45 — 



0.4-4- — 




(a) 



0.12 0.13 0.14- 0.15 0.1S O. 1 7 0.1S O. 1 9 
3 PARAMETER EIT 



0.2 1 




(b) 



0.4-2 0.4-3 0.4-4 0.45 0.46 0.47 0,4-S 0.4-9 0.5 
3 PARAMETER EIT 




(c) 



0.12 0.13 0.14 0.15 0.16 0, 1 7 



Figure 5: The error ellipse projections from which the errors of the 3-parameter fit 
were interpreted, (a) g\ and g2 plane, (b) g2 and plane, and (c) gi and g^ plane. 
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C3 1 



O 



E288 dat^a {PRD 23 (1981) 23} 

CTEQ3. s = (27.4 GeV) = 



5<Q<6 GeV 




g1 = .24. g2= .34, g3= O.. NORM= 0.92 
g1= .15, g2= .48, g3= -.58, NORM= 0.79 



1.2 1.4 

Pt (GeV) 



Figure 6: E288 data, from p + Cu — >■ + X at = 27.4 GeV, with an overall 

systematic normalization error of 25%. The curves are the results of Fit A^i and N2 
and are multiplied by the value of NORM, as shown in the figure and described in 
the text. 

g values are: 
Fit A2 : 

g, = 0.24^°;°? GeV^ = 034;°:°^ GeV^ 
Fit As : 

g, = O.lS^Hs GeV^ , 92 = 0.481°;°^ GeV^ , g^ = -0.581°;^^ GeV'^ ; 

for the 2 and 3 parameter fits, respectively. In Figs. ^ and ^, we show the error ellipse 
projections from which we interpreted the errors of the above fits. In summary, fits 
A2 and ^3 constitute the main results of this paper. 

3.3.3 Cross Checks 

Given these values of g^s and the fitted normalization factors for the E605 data, we 
can calculate the two different predictions for the other three high mass bins not used 
in Fits v42,3, and the results are also shown in Fig. ^. In order to compare with the 
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E288 data, we created Fits A^2,3 in which we fix the g^s to those obtained from Fit 
A3 and fit for NORM from the E288 data alone. Fig. |^ shows the resulting fits are 
acceptable, with values of NORM close to the quoted 25%, namely, NORM= 0.92 
and 0.79 for Fits A'"2,3, respectively. It is encouraging that the quality of the fit for the 
E288 results is very similar to that for the E605 data and that the normalizations are 
now acceptably within the range quoted by the experiment. Hence, we conclude that 
the fitted values of g's reasonably describe the wide-ranging, complete set of data, as 
discussed above. 

We note that although the CDF-Z data, as shown in Fig. |^, contain only 7 data 
points with large statistical uncertainties, they prove to be very useful in determining 
the value of g2- To test this observation, we performed an additional fit, L2. Following 



the method suggested in [0, we set g^ to be zero and fit the gi and (72 parameters 
using the R209 and CDF-Z data alone. (Note that for the R209 data, the typical 
value of T is of the order 0.1, which motivates the choice of tq = 0.1 in the LY form. 
Effectively, the g^ contribution to the R209 data can be ignored.) We found that the 
best fitQ gives g2 = 0.47 GeV^, which is in good agreement with the result of the 
global Fit A3 discussed above. Hence, we conclude that the CDF-Z data already play 
an important role in constraining the g2 parameter, which can be further improved 
with a large Z data samples from Run 1 of the Tevatron collider experiments. We 
shall defer its discussion to the next section. 

4 Run 1 W and Z Boson Data at the Tevatron 

The Run 1 W and Z boson data at the Tevatron can be useful as a test of universality 
and the x dependence of the non-perturbative function W^^{b,Q,QQ,Xi,X2). This 
is clearly demonstrated in Fig. |^, where we give the predictions for the two differ- 
ent global fits (2-parameter and 3-parameter fits) obtained in the previous section 
using the CTEQ3M PDF parameterizations. (The CTEQ4M PDF PB[ gives similar 



results.) With the large Z boson data sets anticipated from Tevatron Run 1 (la and 
5 Also, gi = 0.18 GeV^ 
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Figure 7: Comparison of R209 data with two different theory model predictions 
obtained from the "toy global fit" . The curves are the results of Fit Fi and F2 and 
are multiplied by the value of NORM, as shown in the figure and described in the 
text. 

lb), it should be possible to distinguish these two example models. 

As shown in Ref. [Q, for Qj- > lOGeV the non-perturbative function has little 
effect on the Qt distribution, although in principle it affects the whole Qt range (up 
to Qt ~ <5) • In order to study the resolving power of the full Tevatron Run 1 Z boson 
data in determining the non-perturbative function, we have performed a "toy global 
fit". Fit F3 as follows. First, we generate a set of fake Run 1 Z boson data (assuming 
5,500 reconstructed Z bosons in 24 Qt bins between Qt = Oand 20 GeV/c) using 
the original LY fit resuhs {gi = 0.11 GeV^ ^2 = 0.58 GeV^ and ^3 = -1.5 GeV"^). 
Then, we combine these fake-Z boson data with the R209 and E605 Drell-Yan data 
as discussed above to perform a global fit. The 3-parameter form results in 

Fit F,:g, = 0.10+°;°^ GeV^ , = 0.571^2 GeV^ , 93 = -0.98l°;}^ GeV^^ , 

with a per degree of freedom of approximately 1.3.Q These fitted values for the g^s 

^ This amounts to a shift in the prediction for the mass and the width of the W boson by about 
5MeV and lOMeV, respectively 
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CD 




O 2.5 5 7.5 10 12.5 15 17.5 20 



Qt (GeV) 

Figure 8: Comparison of "fake Z" data (solid curve) with two different theory model 
predictions (dots) obtained from the "toy global fit" Fi and F2. 



Ee05 data (PRD 43 (1991) 2815) 

CTEQ3, s=(38,8 GeV)' 



gl = .24, g2= .34, g3= 0., NORM= 0.94 
gl = .15, g2= .48, g3= -.58, NORM= 0.87 



7<Q<8 GeV 




0.2 0, 



Figure 9: Comparison of E605 data with two different theory model predictions 
obtained from the "toy global fit" . The curves are the results of Fit and F2 and 
are multiplied by the value of NORM, as shown in the figure and described in the 
text. 
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K605 data (PHD 43 (1991) 2815) 

CTEQ3, s = (38.8 GeV)' 



gl= ,24, g2= .34, g3= 0,, NORM= 0.94 
g1 = ,1 5, g?= .48, g3= -.58, NORM= 0,87 




7<Q<8 GeV 



10.5<Q<1 1.5 GeV 



13,5<Q<18 GeV 



Figure 10: Comparison of E288 data with two different theory model predictions 
obtained from the "toy global fit" . The curves are the results of Fit Fi and F2 and 
are multiplied by the value of NORM, as shown in the figure and described in the 
text. 

agree perfectly with those used to generate the fake-Z data except for the value of 
gs, which is smaller by a factor of 2. It is interesting to note that this result agrees 
within 2a with that of Fit ^3, using only the current low energy data, although the 
uncertainties on g^s are smaller by a factor of 2. 

We have also performed the same fit for the 2-parameter form and obtained an 
equally good fit with 

Fit F2:gi = OMt^^m GeV' , g2 = O.eStHa GeV^ 

While the new g2 value (0.63 GeV^) obtained in Fit F2 is very different from that of 
A2 (0.34 GeV^) given in the previous section, it is actually in good agreement with 
the g2 value (0.57 GeV^) obtained from Fit A3. This implies that the Run 1 Z boson 
data, when combined with the low energy Drell-Yan data, can be extremely useful in 
determining the parameter g2. In Figs. |7|-|T0|, we compare the two theory predictions 



(derived from -^2,3) with the R209, fake-Z, E605, and E288 data. As shown, they 
both agree well with all of the data. 

Before closing this section, we would like to comment on the result of a single 
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parameter study of the fake Z data, Fit S\. Given the large sample of Run 1 Z data, 
one can consider fitting the non-perturbative function with only one, Q- independent, 
non-perturbative parameter. With this in mind, we fitted the fake Z data with the 
non-perturbative function 



and found that 



W^;i'^(&,Q,Qo,Xi,X2)=exp -glh' , (15) 



Fit 5i:^ = 2.1+0:09 GeV 



which gives a good description of the Qt distribution of the ''fake-Z" data. It is 
obvious that this fitted value agrees with that of the 2-parameter fit just by considering 
the coefficient of 6^ in first two terms of Eq. |l^. For the results of F2 with the value of 
the Z boson mass, Mz = 91.187 GeV/c^, we obtain 0.04 + 0.63 ln(M^/2Qo) = 2.14, 
which is essentially the same as the coeficient of 6^ in 5*1 . One interesting question 
is whether the result of this one-parameter fit alone can be used to also describe the 
Qt distribution of the boson produced at the Tevatron (at the same energy). 
A quantitative estimate can be easily obtained by noting again that the difference 
between 0.04 + 0.63 ln(Mz/2Qo) = 2.14 and 0.04 + 0.63 ln(MvK/2Qo) = 2.06 with the 
ly-boson mass Mw = 80.3 GeV/c^, is 0.08, which are essentially the same, given the 
uncertainty of 0.09 GeV^ from 5*1. We conclude that it is indeed a good approximation 
to use the one-parameter fit result from fitting Z boson data in order to predict the 
Qt distribution of the boson using the CSS resummation formalism. On the 
other hand, a single parameter without Q dependence (i.e. the parameter ^ alone) 
does not give a reasonable global fit to all of the Drell-Yan data discussed above. For 
instance, for the R209 data, the 2-parameter fit gives 0.04 + 0.63 ln(8/2Qo) = 0.6 for 
the coefficient of Eq. which is not consistent with the value of ^ from Si . Hence, 
we conclude that in order to test the universality of the non-perturbative function of 
the CSS formalism, one must consider its functional form with Q (and x) dependence. 
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5 Conclusions 



The effects of QCD gluon resummation are important in many precision measure- 
ments. In order to make predictions using tlie CSS resummation formalism for the 
Qt distributions of vector bosons at hadron colhders, it is necessary to include contri- 
butions from the phenomenological non-perturbative functions inherent to the formal- 
ism. In this paper, we have extended previous results by making use of 2-parameter 
and the 3-parameter fits to modern, low energy Drell-Yan data. We found that both 
parameterizations result in good fits. In particular our results are 

Fit A2 : 

g, = 0.24;°;°? GeV^ = 0.34l°:°s^ GeV^ 
Fit A3 : 

g, = 0.151°;°^ GeV^ , g, = OASt'of, GeV' , g, = -0.581°;^^ GeV-^ . 

Each functional form predicts measurably different Qt distributions for Z bosons 
produced at the Tevatron. We showed that the full Tevatron Run 1 Z boson data 
can potentially distinguish these two different models. 

In particular, using the results from a toy global fit, we concluded that the large 
sample of the Run 1 Z data can help to determine the value of g2, which is the coef- 
ficient of the ln{Q/2Qo) term in Eq. (0). Given that, one can hope to study the x 
dependence of the non-perturbative function in more detail. We also confirmed that 
it is reasonable to use a single non-perturbative parameter ^ to fit Z boson data, 
and use that result to study the Qt distribution of the boson for Qt < 10 GeV. 
Recently this point has been made in the context of a momentum-space fit [0 us- 
ing a single parameter. Such an approach might indeed alleviate the computational 
overhead required in order to do a complete Fourier transform in order to produce 
distributions of W bosons and decay leptons necessary for analyses. However, if 
one is interested in testing the universality property of the CSS resummation formal- 
ism or making predictions about W and Z boson production at future colliders, such 
as the CERN Large Hadron Collider, then one must include the Q (and, possibly, x) 
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dependent term in the non-perturbative function. 
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